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Abstract 

Several models of dark matter motivate the concept of hidden sectors consisting of 
SU(3)c x SU(2)l x U(l)y singlet fields. The interaction between our and hidden 
matter could be transmitted by new abelian U'(l) gauge bosons A' mixing with 
ordinary photons. If such A"s with the mass in the sub-GeV range exist, they 
would be produced through mixing with photons emitted in decays of rj and rj 
neutral mesons generated by the high energy proton beam in a neutrino target. The 
^4"s would then penetrate the downstream shielding and be observed in a neutrino 
detector via their A' — > e + e~decays. Using bounds from the CHARM neutrino 
experiment at CERN that searched for an excess of e + e~ pairs from heavy neutrino 
decays, the area excluding the 7 — A' mixing range 10~ 7 < e < 10 -4 for the A' 
mass region 1 < Mj^/ < 500 MeV is derived. The obtained results are also used to 
constrain models, where a new gauge boson A interacts with quarks and leptons. 
New upper limits on the branching ratio as small as Br(rj — > 7 A) < 10~ 14 and 
Br{n' — > jX) < 10~ 12 are obtained, which are several orders of magnitude more 
restrictive than the previous bounds from the Crystal Barrel experiment. 
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The understanding of the origin of dark matter has great importance for cos- 
mology and particle physics. Several interesting extensions of the Standard 
Model (SM) dealing with this problem suggest the existence of 'hidden' sec- 
tors consisting of SU(3)c x SU{2)l x U{1)y singlet fields. These sectors of 
particles do not interact with the ordinary matter directly and couple to it by 
gravity and possibly by other weak forces. For example, interaction between 
our and hidden matter may be transmitted by a new abelian U'(l) gauge 
bosons A' (or hidden photons for short) mixing with ordinary photons, first 
discussed by Okun in his paraphoton model [T]. If the mixing strength is very 
weak or the mass scale of a hidden sector is too high, it is experimentally 
unobservable. However, in a class of models the A' may have mass and mixing 
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strength lying in the experimentally accessible and theoretically interesting 
regions. This makes further searches for A n s interesting and attractive, for a 
recent review see [2], and references therein. 

In the Lagrangian describing the photon-hidden photon system the only al- 
lowed connection between the hidden sector and ours is given by the kinetic 
mixing [1112113114] 

L mt = -\<F^ (1) 

where F^, A'^ v are the ordinary and the hidden photon fields, respectively, 
and e is their mixing strength. In the interesting case when A' has a mass, this 
kinetic mixing can be diagonalized resulting in a non-diagonal mass term that 
mixes photons with hidden-sector photons. It means any source of 7's could 
produce kinematically possible massive states A' according to the appropriate 
mixings. Then, when the mass differences are small, photons may oscillate 
into hidden photons, similarly to vacuum neutrino oscillations, with a vacuum 
mixing angle which is precisely e. If the mass differences are large, it could 
results in hidden photon decays, e.g. into e + e~ pairs. 

Experimental bounds on the sub-eV and sub-keV hidden photons can be ob- 
tained from searches for an electromagnetic fifth force [Tf5f6] . from experiments 
using the method of photon regeneration [7l8li9fT0][TT] , and from stellar cooling 
considerations [T21T3] . For example, it has been pointed out that helioscopes 
searching for solar axions are sensitive to the keV part of the solar spectrum of 
hidden photons and the CAST results [l~4"fT5] have been translated into limits 
on the 7 — A' mixing parameter [T6|17|I18|[T9] . Strong bounds on models with 
additional A' particles at a low energy scale could be obtained from astro- 
physical considerations [20]- (22]. However, such astrophysical constraints can 
be relaxed or evaded in some models, see e.g. [23]. New tests on the existence 
of sub-eV hidden photons at new experimental facilities, such, for example, as 
SHIPS [23] or IAXO [25] are in preparation. 

The A n s in the sub-GeV mass range, arising in some models, see e.g. [26|27f28(29 
can be explored through the searches for A' — > e + e~ decays in beam dump 
experiments [30ll31f32ll33f34f35] . or through the rare particle decays, see e.g. 
[36|37f38(39] . For example, if the mass of A' is below the mass of 7r°, it can be 
effectively searched for in the decays n° — > jA', with the subsequent decay of 
A' into e + e~ pair. Recently, stringent constraints on the mixing e in sub-GeV 
mass range have been obtained from a search of this decay mode with existing 
data of neutrino experiments [4"0|4H42] and from SN1987A cooling [43J. 

It should be noted, that many extensions of the SM such as GUTs [33], super- 
symmetric [3S], super-string models [36] and models including a new long- 
range interaction, i.e. the fifth force [31], also predict an extra U (1) factor 
and therefore the existence of a new gauge boson X corresponding to this new 
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group (we denote it X to distinguish from A'). The X's could interact directly 
with quarks and/or leptons, and although the predictions for its mass are not 
very firm it could be light enough (Mx <C Mz) for searches at low energies. 
If the mass Mx is in the sub-GeV range, i.e. of the order of the pion mass, an 
effective search could be conducted for this new vector boson in the radiative 
decays of neutral pseudoscalar mesons P — > 7X, where P = 7r°, rj, or rj', 
because the decay rate of P — > 7 + any new particles with spin or j proves 
to be negligibly small jH] . Therefore, a positive result in the direct search for 
these decay modes could be interpreted unambiguously as the discovery of a 
new light spin 1 particle, in contrast with other experiments searching for light 
weakly interacting particles in rare K, it or /1 decays [4~8|49ll50j . 

Stringent limits on the decay tt° — > 'yA'(X), A'(X) — > e + e~, obtained by using 
results from neutrino experiments have been recently reported in Ref. [3D] , 
The best experimental limits on the branching ratio of the decay 77(7/) — >* lX 
were obtained by the Crystal Barrel Collaboration at CERN [5Tf52] . Using 
proton-antiproton annihilation as a source of X's from 77(77') — > 7X decays, 
they searched for the corresponding 7-peak in their detector, resulting in 90 
% C.L. upper limits on the branching ratio Br(n — > jX) < 6 x 10 -5 for Mx 
masses ranging from 200 to 525 MeV, and Br(i]' — > 7X) < 4 x 10~ 5 for Mx 
between 50 MeV and 925 MeV [5Tf52] . The goal of this Letter is to show that 
more stringent limits on the decay 77(77') — > 'jA'(X), followed by the decay 
A\X) — > e + e~, and the mixing e for A' masses up to ~ 500 MeV can be 
obtained from the results of sensitive searches for an excess of single isolated 
e + e _ pairs from decays of heavy neutrinos in the sub-GeV mass range by the 
CHARM experiment at CERN [531154] . 

The CHARM Collaboration searched for decays — > ue + e~ of heavy neutri- 
nos in the mass range from 10 MeV to 1.8 GeV originated from decays tt, K 
and charmed D mesons decays [53]I54] . The experiment, specifically designed 
to search for neutrino decays in a high-energy neutrino beam, was performed 
by using 400 GeV protons from the CERN Super Proton Synchrotron (SPS) 
with the total number of 2.4 x 10 18 protons on (Cu) target (pot). The CHARM 
decay detector (DD), located at the distance of 480 m from the target, consist 
of decay volume of 3 x 3 x 35 m 3 , three chambers modules located inside 
the volume to detect charged tracks and followed by a calorimeter. The decay 
volume was essentially an empty region to substantially reduced the number 
of ordinary neutrino interactions. The signature of the heavy neutrino decay 
Vf t — > ve + e~ were events originating in the decay region at a small angle with 
respect to the neutrino beam axis with one or two separate electromagnetic 
showers in the calorimeter [SI]. No such events were observed and limits were 
established on the z/ 6jAt — mixing strength as a function of the Vh mass. If 
the decays 77,77' — > ^A' exist, one expects a flux of high energy A n s from the 
SPS target, since neutral mesons 77 and 77' are abundantly produced in the 
forward direction by high energy protons in the target. If A' is a relatively 
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Fig. 1. Schematic illustration of a proton beam dump experiment on search for 
P —7- jA' , A' —7- e + e~ decay chain: neutral mesons P generated by the proton 
beam in the target produce a flux of high energy A' 's through the 7 — A' mixing in 
the decay P — > 77, which penetrate the downstream shielding and decay into e + e~ 
pair in a neutrino detector. The same setup can be used to search for the process 
P -> jX, X -> e+e~. See text. 

long-lived particle, this flux would penetrate the downstream shielding with- 
out significant attenuation and would be observed in the CHARM detector via 
the A' — > e + e~ decay into a high energy e + e~ pair, as schematically illustrated 
in Fig. [H The occurrence of A' — > e + e~ decays would appear as an excess of 
e + e~ pairs in the CHARM DD above those expected from standard neutrino 
interactions. The experimental signature of these events is clean and they can 
be selected in the CHARM DD with a small background. As the final states 
of the decays v h — > z/e + e~ and A' — > e + e~ are identical, the results of the 
searches for the former can be used to constrain the later for the same e + e~ 
invariant mass regions. 

The flux of hidden photons from decays of r/'s and r/'s produced in the target 
by primary protons can be expressed as follows: 

, ,, f d 3 a(p + N — > 77(77') + X) , , .. o 
Q(A') ex N pot j KP * 1 ] ^Br^r/) n)f#Mtf) ( 2 ) 

where N pot is the number of pot, a(p + iV — > r](j]') + X) is the f](r]') me- 
son production cross-section, Br(r](r]') — > 77) is the r](r]') — > 27 decay mode 
branching fraction [55], and / is the decay phase space factor, respectively. 

To perform calculations we used simulations of the process shown in Figfj] 
from our previous work on tt° — > 'jA' decays [40J by taken into account the 
relative normalization of the yield of different meson species 7r° : 77 : vj from 
the original publications. The invariant cross section of a hadron production 
can be expressed as [55] 

d 3 p pTdprdydcj) 1-npTdpTdy 

where px is the transverse momentum of the particle, y is its rapidity, and in 
the last equality integration over the full 2tt azimuthal angle <fi is performed. 
For the production cross sections of 77 and 77' neutral mesons we used the 
Bourquin-Gaillard (B-G) formula from Ref.[56], which gives the parametric 
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Fig. 2. Invariant cross section for n production in pp collisions at y/s = 27.4 GeV 
vs pt and rapidity y c . m . over the range —3.0 < y c .m. < 3.0. 

form of ([3]) for the production in high-energy hadronic collisions of many 
different hadrons over the full phase-space: 

d 3 a(p + N ^7](r]')+X) . ( 2 xi2.s 5.13 . , , AS 
L = A ^(e^2) ex P(-^ )/( ^ } ' (4) 

where 

exp(-pT) Pt < 1 GeV/c 

exp(-l - 23(p T - l)/y/s) p T > 1 GeV/c 

with Y = y m ax — y, being the rapidity, and Et(pt) the transverse energy 
(momentum) in GeV. Coefficients A^n in (4) are normalization factors that 
were tuned to obtain the cross sections a v and a v > in pp collisions at 400 
GeV/c. In these calculations for the relative yields of 7r° : rj : rf mesons the 
values 1 : 0.078 : 0.024, respectively, were used, which were obtained from 
the measurements of the total cross sections a n o (127.2 ± 1.5 ± 3.2 mb) and a v 
(9.78 ±0.56 mb), and from an estimate of the rf production in pp-interactions 
at 400 GeV/c by the NA27 experiment at CERN SPS [57]. The invariant rj 
production cross section, obtained by using the B-G parameterization (4), is 
shown in Fig. [2] as a function of pt and rapidity y c . m . in the center of mass 
system. The total r] and rf production cross sections in p-Cu collisions were 
calculated from its linear extrapolation to the target atomic number. The B- 
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G px distributions at low transverse momenta were corrected by taking into 
account the measurements results of 7r°'s and 77's production in 450 GeV p- 
Be and p-Au collisions down to very low p T (> 20 MeV/c) obtained by a 
joint experiment of the TAPS and CERES Collaborations at the CERN SPS 
[58] . In this experiment precise measurement data on //-production at low pt, 
revealed a lower cross section than expected from extrapolating data at higher 
Pt- The B-G parameterization, found to be in a good agreement with tt° data 
for the full pt range, and also with the rj data for the large (pt > 0.4 GeV/c) 
momentum transfer, overestimates the production rate of rj mesons by 25-30% 
at transverse momenta pt < 0.4 GeV/c. This observation results in decrease 
of the acceptance of the CHARM DD for A n s produced in 77(7/) — > ^A' decays 
and, hence has to be taken into account. To achieve a better description of 
the \cw-pt -region, a parameterization suggested by [581159] was used: 

where ttit is the transverse mass, xt = 2mr/ y/s, and parameters — 0.15, b = 
6.5 GeV -1 , a = 0.011, c = 7.9 and normalization factors were deter- 

mined in Ref. [58] by fitting the precise low pt experimental data in the rapidity 
range 3 < y < 4. The results of fitting were also in good agreement with other 
meson production data obtained for different rapidity intervals. The compar- 
ison of two parameterizations is shown in Fig. [3j We used Eq.((5]) also for 
calculations of low pt distributions of 7/ mesons. 
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Fig. 3. Invariant cross section for 77 mesons production in pp collisions at yfs = 27 A 
GeV calculated as a function ofpT in the Iow-pt region for the B-G parameterization 
of Eq.(4) (solid curve), and from Eq.^( dashed curve). 
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The calculated fluxes and energy distributions of r], rj' produced in the Cu 
target were used to predict the flux of A"s, as a function of its mass. For a 
given flux cL$(Ma', Ea>, N pot )/dEA' of A n s the expected number of A' — > e + e~ 
decays occurring within the fiducial length L of the CHARM detector located 
at a distance V from the neutrino target is given by 



N 



/d<§> 

uEa 1 

C,AdE A > 



•exp 



L'Mj 
Pa'T A > 



1 — cxp 



LM A 
Pa'Ta 



(6) 



where Ea 1 , Pa 1 -, and ta> are the A ! energy, momentum and the lifetime at rest, 
respectively, and ( is the e + e~ pair reconstruction efficiency. The acceptance A 
of the DD was calculated tracing A n s produced in the Cu-target to the detector 
taking the relevant momentum and angular distributions into account. As an 
example for a mass Ma> = 300 MeV, A = 0.07 and ( ~ 0.6 [53J. In this 
estimate the average momentum is < pa> >— 35 GeV and L ~ 10 m. 

The obtained results can be used to impose constraints on the previously 
discussed models with A' hidden photons. For A' masses below the mass M^^ 
of 77(77') meson, the corresponding branching fraction of the decay 77(77') — > 7 A', 
is given by [26|l27] : 



Br(r] -)• 7^4') = 2e 2 Br(r] ->• 77) ( 1 - 



Ml ^ 



M 2 



(7) 



with the similar expression for the 77' decay. Assuming the dominant A' boson 
decay is into e + e~ results in a corresponding decay rate, which for small mixing 
is given by: 



T{A' ->• e+e" 



Q 



e 2 M A < 



Ami ( 2m 
1-^(1 + 



M\, 



Using the relation N®+f_ > N A >-> e +e-, where N^L (= 2.3 events) is the 
90% C.L. upper limit for the number of signal events pj3Tol] and Eqs.(7,8), 
one can determine the 90% C.L. exclusion region in the (ikf^e) plane from 
the results of the CHARM experiment, which is shown in Fig. H] together with 
regions excluded by the Nomad [10] and electron beam dump experiments 
E137, E141, E774 [3lpll33lM] . and by recent measurements from APEX [35], 
KLOE [361 , BaBar EH1 and MAMI (391. 



The shape of the exclusion contour from the CHARM experiment correspond- 
ing to the A' mass range Ma> ^ 300 MeV is defined mainly by the phase 
space factor in (JTj). Using similar considerations the exclusion area for the 
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Fig. 4. Exclusion region in the (M^r,e) plane obtained in the present work from 
the results of the CHARM experiment |53|5^y. The areas excluded from the (g-2) 
considerations, the results of the Nomad fjm and electron beam dump experiments 
E137 JSfflSiy . E141 JMI, E774 and from the searches by APEX JM$, KLOEJM/, 
BaBarlW^ and MAMI are also shown for comparison. 

decay 7/ — > jA' has been also obtained. However, due to the lower 7/ produc- 
tion cross section and a small branching fraction of the decay 7/ — > 77 this 
exclusion area falls basically within that of obtained for rj decays. Note, that 
the uncertainty for the rj production at low does not significantly affect the 
limits derived. For example, the variation of the r\ yield in (j6]) by 25% results 
in the corresponding variation of the limits on mixing strength e of the order 
of 5%. This is because the sensitivity of the search is proportional to the e 4 . 
Indeed, in ([6]) the branching fraction of ([7]) and the decay rate T(A' — > e + e~) 
of ([8]) both are proportional to e 2 . 

We can also constrain models where the previously discussed X bosons interact 
with both quarks and leptons. The 90% C.L. upper limits on the Br(r](r]') — > 
r )X)Br(X — > e + e~) vs X lifetime tx shown in Fig. [5] were calculated by 
using the relation > N x ^ e +e~ and Eq.([6]). Our result is sensitive to a 
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branching ratio Br(f](i]') — > 7X) > 10~ 14 (10~ 12 ), which is about nine (eight) 
orders of magnitude smaller than the previous limit from the Crystal Barrel 
experiment [51 |I52"] . Over most of the t x region, the X lifetime is sufficiently 
long, that LM X /px T x *C L'M x /px T x "C 1. The interaction of X bosons with 
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Fig. 5. The 90% C.L. upper limits on the branching ratio 
Br(j]{rf) — > jX)Br(X — > e + e~) versus tx obtained from the CHARM ex- 
periment for rj (solid curves) and rf (dashed curves) decays. The numbers near the 
curves indicate the corresponding values of Mx ■ 

both quarks and leptons can be written in the form: 

L x = 9x{QbxB 1 + Q eX Ll + Q^xLl + Q rX L l T )X l (9) 

where «j = gjc/^ 71 is the coupling constant, B l = S ? = u ,d,s,.. qYQi L l e = eYe + 
VeLY u eL, see e.g. [I81I50] . Assuming charges Qbx — Q e x — 1, we found 

1 / Ml \ - 3 / 2 
ax< k (l-f!^) (10) 

M x [MeV]\ M^J 1 ; 

which is valid for k = 4.3 x 10" 12 and M x < 2M W ( 270 MeV). For the mass 
range from 270 to M v (548 MeV) the decay channel X — > tttt is open. To 
calculate this decay rate is not a simple task, because the QCD long distance 
effects have to be taken into account [IS]. Therefore, to avoid this difficulty, 
we used for this decay rate a reasonable estimate T(X — > nn) ~ cyxQ 2 bx^x, 
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which results in k = 1.4 x 10 -12 for Qbx — 1- From the similar consid- 
erations, for the mass range M v < Mx < M r? /(948 MeV) the limit on the 
coupling constant is given by fllOp with replacements k = 3.5 x 10~ 10 and 
Mrf instead of M v . The bounds ([TO]) are valid for r x > 10~ 10 M x [MeV] s. 
They are more restrictive than those obtained in [49|50] , and than bounds re- 
ported by NOMAD [12]. Less stringent limits (by a factor ~ a) could be 
obtained for the cases where the X interacts only with leptons, or when 
it is a leptophobic boson which interacts only with quarks and decays into 
e + e~ pair through the quark loop [4"9"f5D] . For X produced in r\ decay with 
the masses M x < 545 MeV/c 2 the best limits from CHARM are in the 
region Br(i] — > r yX)Br(X — > e + e~) < (2 — 3) x 10~ 14 , while the corre- 
sponding limits for X's with masses Mx < 948 MeV/c 2 from rf decays are 
Br(rf -> "fX)Br(X -> e + e~) < 2 x 10~ 12 . The attenuation of the X- flux due 
to X interactions with matter was found to be negligible, e.g. for couplings of 
(11 01) the X boson mean free path in iron is much longer, as compared with the 
iron and earth shielding total length of 0.4 km used for the CHARM beam. 

In summary, using sensitive limits from the CHARM experiments on heavy 
neutrino decays v h — > ue + e~ , new bounds on a hidden sector gauge A' bosons 
produced in the neutral meson decay f](r]') — > 7 A' in the sub-GeV A' mass 
range are obtained. The A n s could mediate interaction between our world and 
hidden sectors consisting of SU(3)c x SU(2)l x U{1)y singlet fields through 
the mixing with ordinary photons. For the A' mass region 1 < Ma> ^ 500 
MeV the obtained exclusion area covers the 7 — A' mixing strength in the 
range 10~ 7 < e < 10~ 4 . The results obtained are also applicable to the rj 
and rj' decays into a photon and a new gauge bosons X that interacts with 
quarks and leptons, or only with quarks. Our best result is sensitive to a 
branching ratio Br(r](ri') — > 7X) > 10~ 14 (10~ 12 ), which is about nine (eight) 
orders of magnitude stronger than the previous limit from the Crystal Ball 
experiment [5 1 |I52"] . The obtained constraints were used to set new limits on 
the X coupling strength to lepton and quarks which are more restrictive than 
bounds from previous searches. These results enhance existing motivations for 
further more sensitive search for A'(X) decay at the high intensity frontier 
[60] and additional tests of hidden sectors in neutrino experiments. 

The help of V.D. Samoylenko and D. Sillou in calculations is greatly appreci- 
ated. 
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